Magnetic structures produced by the small-scale dynamo 
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Small-scale dynamo action has been obtained for a flow previously used to model fluid turbu- 
lence, where the sensitivity of the magnetic field parameters to the kinetic energy spectrum can be 
explored. We apply quantitative morphology diagnostics, based on the Minkowski functionals, to 
magnetic fields produced by the kinematic small-scale dynamo to show that magnetic structures are 
predominantly filamentary rather than sheet-like. Our results suggest that the thickness, width and 
length of the structures scale differently with magnetic Reynolds number as i^m 2 ^ 1- ^ and R^ ' 55 
for the former two, whereas the latter is independent of i? m , with s the slope of the energy spectrum. 



o 
o 

<D 
CO 

<3\ 



<N 
> 

<N 
<N 
O 
t> 

o 

Oh! 

o 



PACS numbers: 47.65.-d, 07.55.Db, 02.40.Pc 

The fluctuation (or small-scale) dynamo is a turbulent 
dynamo mechanism where a random flow of electrically 
conducting fluid generates a random magnetic field with 
zero mean value (the other type of dynamo relies on de- 
viations of the flow from mirror symmetry, and perhaps 
other symmetries, and produces mean magnetic fields). 
This dynamo mechanism appears to be responsible for 
the random magnetic fields in the interstellar medium 
[H and in galaxy clusters 0, 0, Q|. A necessary con- 
dition for the small-scale dynamo is that the magnetic 
Reynolds number R m is large enough that the random 
velocity shear dominates over the effects of the fluid's 
electric resistivity. Here R m = ul/rj where u and I are 
a typical velocity and length-scale respectively and n is 
the magnetic diffusivity. The critical magnetic Reynolds 
number, Rip er, is 30-500, depending on the nature of the 
flow 0, @, 0, 0, [Io|, [H| • Numerical simulations of the non- 
linear, saturated states of the small-scale dynamo have 
been reviewed in ref. pjj , but more insight can be gained 
by a careful analysis of the kinematic regime where the 
magnetic field is still too weak to affect the flow. In par- 
ticular, the turbulent inert ial range is undeveloped for 
the modest values of the Reynolds number achievable in 
direct numerical simulations; therefore, numerical results 
are often ambiguous. 

Here we study the kinematic small-scale dynamo by 
solving numerically the induction equation 



-=Vx(uxB) 



r?V 2 B, V-B = 0, (1) 



for the magnetic field B in a prescribed flow u(x, t). Un- 
like other numerical models of the small-scale dynamo, 
our choice of the velocity field allows us to control fully 
its energy spectrum and time variation. In particular, 
the flow has a well pronounced power-law spectral range 
with controllable spectral slope. We use a spectral model 
for u, specified in Eq. ([2]), which was developed as a La- 
grangian model of turbulence and contains its essential 
features: it is three dimensional, time-dependent, multi- 
scale and has transport properties (two-particle disper- 
sion, etc.) which agree remarkably well with those of tur- 



bulent flows HI, j3] in both experiments and numerical 
calculations. Therefore, our main results are comparable 
with numerical experiment s (lR \v\ and to some extent 
with experimental results [15]. The velocity field u is 
given by 
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where <\> n = k n -x+cc; n t, N is the number of modes and k n 
are randomly chosen unit vectors (k n = fc n k n is the wave 
vector); note that u is solenoidal by construction. The di- 
rections of C n and D n are chosen randomly; we require, 
however, that they are normal to k n , so that the root 
mean velocity of each mode is [{C 2 +D 2 l )/2} 1 / 2 . The mag- 
nitudes of C n and D n are chosen to reproduce the desired 
energy spectrum, E(k): C n = D n = [^E(k n )Ak n } 1/2 , 
where Ak n = ^(fc n +i — fc n _i) for 2 < n < N — 1, but 
Aki = |(&2 — fci) and Afcjv = \{kw — fcjv-i)- For 
example, a model spectrum with suitable behaviors at 
k — > 0, k — > 00, and E(k) ex k s within the inertial range 
fco <C k <C fed, can be implemented with 
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where fco is the integral scale of the flow, fcd is the dissi- 
pation wavenumber and a is the normalization constant 
used to control the intensity of the flow. The Kolmogorov 
inertial range spectrum is obtained for s = —5/3. The 
frequencies co n = [fc^(fcn)] 1 / 2 introduce time variation 
such that each mode varies at its 'eddy turnover' time. 

Having specified fco and fcd, we can introduce the ef- 
fective Reynolds number via Re = (fcd/^o)^ 1_s ^ 2 and 
R m = (u 2 ) 1 ^ 2 lo/r] for the magnetic Reynolds number, 
where angular brackets denote suitable averaging and 
lo = 27r/fco- We enforce periodicity in the unit box by 
choosing k n randomly from a family of vectors whose 
components are multiples of 2tt. The resulting kinetic 
energy spectrum E{k) for s = —5/3 is shown in Fig. [T] 
(solid line), defined as f™E(k)dk = V' 1 f v ±\u\ 2 dV, 
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FIG. 1: Typical energy spectra for the velocity field (E, solid) 
and a growing magnetic field (M, dotted) in arbitrary units, 
with N = 40, ho = 2tt and kd ~ 167r; this corresponds to the 
effective Reynolds number Re « 16 and Rm ~ 2000. Dashed 
lines represent power laws with the slopes indicated. 
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FIG. 2: The growth (or decay) rate of the r.m.s. magnetic 
field versus the magnetic Reynolds number for a flow shown 
in Fig.[T]with Re w 16; note that a = for R m — R m ,cr ~ 753. 
The dashed curve is an asymptotic fit as discussed in the text. 



where V is total volume, as obtained numerically via 
Fourier transform of u given by Eq. ([2j) with the input 
from Eq. (j3j. Although the flow is not random, it exhibits 
Lagrangian chaos 



13L Il4j , which facilitates the dynamo 



action. The flow is stationary on average, so we expect 
that the mean magnetic energy density will vary expo- 
nentially in time, (B 2 ) 1 / 2 oc e at , with the growth (or 
decay) rate a depending on R m and the kinetic energy 
spectrum [7]. Meanwhile, the spatial form of the mag- 
netic field is expected to remain statistically stationary. 

We solve the induction equation (pQ) in a unit cubic pe- 
riodic box discretized in 128 3 grid cubes. We begin with 
a periodic, seed magnetic field with energy distributed 
equally at 5 different scales in the box. After a small 
number of timesteps, the 5-peaked magnetic energy spec- 
trum smooths out to become a horizontal fiat spectrum. 
Figure [2] shows a, measured in terms of the turnover time 
at the smallest velocity scale, 27r/cu n . Dynamo action 
is initiated at i? mjC r ~ 753 for the Kolmogorov spec- 
trum, s — —5/3. The data were obtained by chang- 
ing R m (or rj) without changing the flow. Thus, the ef- 
fective magnetic Prandtl number P m = R m /Re varies 
together with R m . Also plotted (dashed curve) is the 



asymptotic behavior suggested for R m ^> 1 in ref. 0]: 
a « a(U/C) In (1Z m /1Z m ^ cr ), where U and C are a repre- 
sentative speed and length respectively and a is a con- 
stant of order unity; [7] obtain a ~ 2/3 for U = uq, 
C = 27r//co- A very good fit to our results is achieved for 
a similar value a ~ 0.6 (or w 0.2 if the dissipation scale 
is used). The accuracy of this fit for R m > R m , C r is over 
96%. Magnetic structures at the largest values of R m and 
k considered here may be only marginally resolved with 
the resolution 128 3 (cf. ref. [12]) but note that the size of 
our computational box is unity rather than 2ir. However, 
the systematic behavior of the growth rate and magnetic 
spectrum at the largest values of R m and k indicate that 
the numerical resolution is sufficient; we carefully checked 
for other signs of insufficient resolution in our results but 
none were found. 

The energy spectrum M(k) of the growing magnetic 
field (defined similarly to E(k)) for R m w 2000 is shown 
in Fig. [1] (dotted line), where it exhibits a range with 
M(k) oc k^l 2 expected for a small-scale dynamo with a 
single-scale random flow [5]. As we discuss below, our 
model also reproduces all other known features of the 
small-scale dynamo that we have tested. 

Asymptotic solutions of the induction equation for 
R m ^> 1 relevant to the small-scale dynamo have been 
interpreted as indicating that magnetic field is mostly 



concentrated into filaments of thickness l v = IqRj 



-1/2 



(for single-scale flows) [6] or = loRm 2 ^ 1- ^ for flows 
with kinetic energy spectral index s @, H, 0]. However, 
more recent numerical simulations appeared to display 
magnetic sheets and ribbons To|, EH E3 . The latter con- 
clusion was based on the visual inspection of magnetic 
isosurfaces and application of heuristic morphology indi- 
cators. Here we apply mathematically justified morphol- 
ogy quantifiers, based on the Minkowski functionals [lil ]. 
to isosurfaces of magnetic energy density such as those 
shown in Fig. 03 This tool has previously been applied to 
galaxy distribution and cosmological structure formation 

The morphology of structures in three dimensions can 
be fully quantified using the four Minkowski functionals: 



V 2 = // («i + k 2 ) dS , V 3 = j j kik 2 dS , 



(4) 



where integration is over the volume and surface of the 
structures, respectively, and n\ and K2 are the prin- 
cipal curvatures of the surface. The Minkowski func- 
tionals can be used to calculate the typical thickness, 
width and length of the structures, as T = Vq/2Vi, 
W = 2V 1 /7rV 2 , and L = 3V 2 /4V 3 , respectively. Then, 
useful dimensionless measures of 'planarity' P, and 'fila- 
mentarity' F [19] can be denned as P = (W-T)/(W+T), 
F = (L — W) /{L + W). In idealized cases and for convex 




FIG. 3: (color online) The B 2 = 3{B 2 ) isosurfaces of magnetic 
energy density for the flow of Fig. [T]with R m = 1570. 



surfaces, values of P and F lie between zero and unity. 
For example, an infinitely thin pancake has (P, F) = 
(1,0), a perfect filament has (P,F) = (0,1), whereas 
(P, F) = (0, 0) for a sphere. A B 2 = 1.5(P 2 ) isosurface of 
the initial magnetic field has (P, F) = (0.094, 0.14). Some 
other examples are shown on the right-hand side of Fig.[H 
We note that the unit cube has T = 3/4, W = 2/tt, 
L = 1/2, thus these measures are not necessarily such 
that T < W < L. Deviations from this ordering are rel- 
atively rare for random fields studied here, yet to avoid 
confusion we introduce the notation l\ = min(T, W,L), 
l 2 = med(T, W, L) and l 3 = max(T, W, L). 

As shown in Fig.JU F increases faster than P with P m , 
so the magnetic structures produced by the small-scale 
dynamo are better described as filaments, especially at 
the larger values of R m . Remarkably, similar velocity 
field structures are not filamentary since the wave vec- 
tors in Eq. ([2]) have no preferred direction, hence the 
velocity field at small scales is nearly isotropic. Cor- 
respondingly, the isosurfaces of u 2 have negligible pla- 
narity and filament arity. The isosurface of vorticity ft 
with ft 2 = 4(ft 2 ) has (P,P) = (0.18,0.11); similarly the 
isosurface of the total strain (S 2 = SijSij) with S 2 = 
4.5(S 2 ) has (P,P) = (0.11,0.16). Thus, the morphology 
of the magnetic field is controlled by the nature of the 
dynamo action rather than by immediate features of the 
velocity field. Isosurfaces of the electric current density 
J = V x B are ribbon-like, with (P, F) = (0.57, 0.82) at 
a level J 2 = 4( J 2 ) 2 for R m = 1570. 

Using the Minkowski functionals, we can also reli- 
ably measure the characteristic length scales of magnetic 
structures and explore their scalings with R m and s. In 
Fig. [5j we display h against R m at two instants in time 
for the flow with s = —5/3 shown in Fig[TJ Whilst the 
behavior for R m <C P m ,cr shows variations in time, we 
observe for R m > 200 a time-independent scaling of the 
thickness of magnetic structures: l\ = 2ir/k ri ~ Rm 3 ^- 
This scaling, obtained in ref. for a flow ^-correlated in 
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FIG. 4: (color online) The route taken by magnetic structures 
through the (P, P)-plane as Rm increases, averaged for each 
R m between eight isosurfaces ranging from B 2 = 1.5(P 2 ) to 
B 2 = 5(P 2 ) in steps of 0.5(P 2 ); in this range the magnetic 
structures depend only weakly on the choice of isosurface. 
Three data points are labelled with their values of R m - The 
dashed line is the estimates of Eq. J6]). The objects on the 
right, from a filament to a sphere, are shown for reference: 
= (a) (0.096,0.81); (b) (0.66,0.23); (c) (0.66,0.12); (d) 
(0.25,0.66); (e) (0.18,0.43); (f) (0.14,0.23); (g) (0.087,0.073) 
and (h) (0.0036, —0.0047). The cross section (in green) is also 
displayed for each object, except for (h) where it is circular. 
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FIG. 5: Average values of the thickness h of the isosurface 
structures against R m at two instants in time for the flow 
of Fig. [T] Data marked 'x' ('+') are results attained when 
the smallest eddies have made 25 (50) revolutions. The solid 
straight line has a gradient of —3/4. 

time (see also ref. 0, [s|) follows from the balance of the 
induction and diffusion terms in the induction equation 
[22I ]: a flow with energy spectrum E(k) oc k s yields 

Zi ^ 0J R- 2/(1 - s) . (5) 

The reason for the i? m - dependence of this scale is akin to 
that of the magnetic diffusion length, or the skin-depth. 

Figure [6] shows the characteristic width I2 and length 
Is of the magnetic structures against R m . For R m > 
200 we observe another time-independent scaling, I2 ~ 
i?~ 0,55 . This distinct behavior of the width of magnetic 
structures has not been obtained in earlier analytical or 




4 



0.1 


1 1 — 1 — 1 — 1 — 1 — 1 1 1 1 — 


~l — 1 — 1 — 1 — 1 — 1 1 

X 








fx 

frX 
I X 

Vx 


n ni 

U.U I 


0.1 

r ' 3 

0.01 


E i i ■ E 




I 


100 1000 

1 1 1 


1 



100 „ 1000 



FIG. 6: As in Fig. [5] but for the width l 2 (main frame) and 
length /3 (inset). The solid straight line in the main plot has 
a gradient of —0.55. 

numerical studies of the small-scale dynamo. The simul- 
taneous decrease of I2 and h : coupled with the approxi- 
mately i? m -independent behavior of £3 (Fig. [6] inset) sup- 
ports the notion that the magnetic structures become 
filamentary as R m increases (see Fig. |4]). Indeed, using 
h ~ 2.4i?-°- 75 , l 2 ~ 0.9i?-°- 55 and l 3 ~ 0.05^ in the 
definitions of P and F, we can estimate, for s = —5/3, 

P-1-2I1 + K 2 ]- 1 , F~l-2[l + ^< 55 ]" 1 , (6) 

so that F > P for R m > 200. These relations are shown 
by the dashed line in Fig. [H 

To investigate how the scaling laws identified via the 
Minkowski functionals compare with those inferred from 
other measures, we calculated the inverse 'integral scale' 
of the magnetic field, 2tt/Zj = / kM{k)dk/ J M(k)dk. 
We found that lj follows a scaling of ^~ 0,42 , which un- 
derstandably differs from the scalings of /1, I2 and £3. The 
scale lj is a poor measure of the dimension of anisotropic 
magnetic structures such as filaments. We note that the 
above scaling of lj is maintained for all subcritical and 
supercritical values of i? m , unlike the results of Figs. [5] 
and [6] which display well-defined, time-independent scal- 
ings only for R m > 200. We have verified that the scal- 
ing ([5]) emerges for s = —5/3, —2, —3. On the contrary, 
I2 ~ ^m°' 55 independently of s. Asymptotic solutions 
[3] suggest that the small-scale dynamo (with Pr m < 1) 
is only possible for s < —3/2. Our results show that, 
for high effective Pr m , the dynamo action is possible for 
s = — 1 as well, although a scaling different from ([5]) is 
exhibited. In this context, we propose that the nature of 
the asymptotic solution, rather than the possibility of a 
dynamo, is different at s = — 1 from that at s < —3/2. 
We also tested the proxy dimensions of ref. [23] but these 
do not show physically justifiable scalings. 

The scaling of l\ with R m has been obtained @, H, 
for flows with Pr m = R m /Re < 1, where the maximum 
of the magnetic spectrum occurs in the inertial range 
of turbulence. With the effective Re defined above, our 



simulations have an effective Pr m > 1. We cannot intro- 
duce a true Re because we do not solve the Navier- Stokes 
equation and therefore the kinematic viscosity is not de- 
fined in our model. The reason why scalings expected 
for Pr m < 1 occur in our model which appears to have 
Pr m > 1 remains unclear. 

In conclusion, we have shown that velocity fields of the 
form ([2]) are able to initiate the small-scale dynamo. This 
model flow reproduces many properties of turbulence and 
thus offers a convenient tool to study dynamos in multi- 
scale, well controlled flows. A quantitative, mathemat- 
ically justifiable description of the morphology of mag- 
netic structures based on the Minkowski functionals in- 
dicates that, at least at the kinematic stage, magnetic 
field is concentrated into filaments rather than sheets 
or ribbons. We have confirmed that the characteristic 
thickness of magnetic structures scales with the mag- 
netic Reynolds number as given in Eq. ([5]) for the energy 
spectrum slope of s, whereas the width varies as R^ ' 55 
independently of s, and the length is independent of R m . 
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